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We derive and analyze the Bethe-Salpeter equations for spin zero diquarks in the color super- 
conducting phase of cold dense QCD with two massless flavors. The spectrum of diquarks contains 
an infinite number of massive excitations and five (nearly) massless pseudoscalars. The former are 
singlets while the latter include a doublet, an antidoublet and a singlet with respect to the unbroken 
^^ , SU{2)c- Because of an approximate parity doubling at large chemical potential, all massive states 

come in pairs. The decay constants, as well as the velocities of the (nearly) massless pseudoscalars 
are derived. The different role of the Meissner effect for tightly bound states and quasiclassical 
bound states is revealed. 
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""' I. INTRODUCTION 

^ , With continuing advances in modern nuclear and high energy experiments, it has become feasible to produce 

^^ ' deconfined quark matter in the laboratory. Not surprisingly, this has stimulated many theoretical studies of quark 
^_" , matter at high densities and/or temperatures (for recent reviews see, for example, Refs. [|l|-Ql). 

r^ ' Of special interest is the cold quark matter at high densities, i.e., at densities which are at least a few times 

f**) ' larger than that of a nucleon no ~ 0.17 fm~^. It was known for a long time that such matter could be a color 

f^ t superconductor pUq]. Nevertheless, until recently our understanding of the color superconducting state had remained 

^^ ' very poor. The new developments started with the ground breaking estimates of the color superconducting order 

(-H , parameter in Refs. nn- Within a phenomenological (instanton liquid) model, it was shown there that the order 

p) ' parameter could be as large as 100 MeV. These estimates revived the hope of producing and detecting the color 

' , superconducting phase either in experiments or in natural systems such as neutron (quark) stars. Being motivated 

TT^' by the potential possibility of observing the color superconducting phase at moderate densities, the subject bursted 

(-H [ with numerous studies and new discoveries. 

• • . At first, it was not clear at all that the microscopic theory, quantum chromodynamics (QCD), would lead to the same 

. !^ ' (or, at least, the same order of magnitude) estimates for the superconducting gap as the phenomenological models. It 

^ . was suggested in Refs. |^,|l^ that the screening effects of gluons should play a crucial role in the analysis. In particular, 

^ ' while the electric gluon modes are subject to the Debye screening at already the scales of order Id ^ i/gslL, where fj. 

. 5t 1 is the chemical potential and gs is the running coupling related to the scale /i, the magnetic modes are subject only 

to the Landau damping, which does not completely eliminate the long range interaction ||lO(| . The subsequent studies 

of the gap equation in QCD confirmed that proper treatment of the gluon screening effects are crucial in deriving the 

estimates for the superconducting gap ||ll|-[l7|] . 

Also, it was revealed that the ground state of the quark matter with three light flavors is given by the so called Color 
Flavor Locked (CFL) phase ||l^. It is remarkable that the chiral symmetry in such a phase is spontaneously broken 
and most of the quantum numbers of physical states coincide with those in the hadronic phase. It was tempting, 
therefore, to suggest that there might exist some kind of continuity between the two phases [^. Another interesting 
feature of the three flavor QCD was pointed out in Ref. [Q, where the possibility of gapless color superconductivity (a 
metastable phase) was proposed. In addition, many interesting patterns of symmetry breaking were revealed in models 
with the number of flavors larger than three [[l7| , as well as in two-color QCD with the quarks in the fundamental 
representation and in any-color QCD with quarks in the adjoint representation pit]. The anomaly matching conditions 
were analyzed in Ref. l22| . 
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The low energy dynamics of the color superconducting phase could be efhciently studied by using effective actions 
whose general structure is fixed by symmetries ||23H26|]. The finite set of parameters in such theories could be 
either taken from an experiment (when available), or sometimes derived from QCD (for example, in the limit of 
the asymptotically large chemical potential) . Because of the nature of such an approach, at best it could probe the 
properties of the pseudo-NG bosons, but not the detailed spectrum of the diquark bound states (mesons). It was argued 
in Ref. [g^j however, that, because of long-range interactions mediated by the gluons of the magnetic type P,[lO|, the 
presence of an infinite tower of massive diquark states could be the key signature of the color superconducting phase 
of dense quark matter. 

In this paper, we consider the problem of spin zero bound states in the two flavor color superconductor using the 
Bethe-Salpeter (BS) equations. (A brief outline of our results was given in Ref. |Q). We find that the spectrum 
contains five (nearly) massless states and an infinite tower of massive singlets with respect to the unbroken SU{2)c 
subgroup. Furthermore, the following mass formula is derived for the singlets: 



M,^-4|A,^l%l-(^;^j, n=l,2,..., (1) 

where k is a constant of order 1 (we find that k ~ 0.27), | Aq | is the dynamical Majorana mass of quarks in the color 
superconducting phase, and as = gl/^'K. 

The Meissner effect plays a crucial role in obtaining this result. In particular, the important point is that while the 
Meissner effect is essentially irrelevant for tightly bound states, it is crucial for the dynamics of quasiclassical bound 
states (whose binding energy is small). 

At large chemical potential, we also notice an approximate degeneracy between scalar and pseudoscalar channels. 
As a result of this parity doubling, the massive diquark states come in pairs. In addition, there also exist five massless 
scalars and five (nearly) massless pseudoscalars [a doublet, an antidoublet and a singlet under SU{2)c\- While the 
scalars are removed from the spectrum of physical particles by the Higgs mechanism, the pseudoscalars remain in the 
spectrum, and they are the relevant degrees of freedom of the infrared dynamics. At high density, the massive and 
(nearly) massless states are narrow resonances. 

This paper is organized as fol low s. In Sec. O, we describe the model and introduce the notation. Then, further 



developing our notation in Sec. Ill, we briefly review the approach of the Schwinger-Dyson equation in the color 
superconducting phase of Nf — 2 QCD. In Sec. N^, we derive the Ward identities for the quark-gluon vertex functions, 
corresponding to the broken generators of the color symmetry. These identities are going to be very helpful in the rest 
of the paper. We outline the general derivation of the Bethe-Salpeter equations for the diquark states in Sec. M. The 
detailed analysis of the Bethe-Salpeter equations for the NG bosons and the massive diquarks is presented in Sees. 



and VII, respectively. Appendix R^ contains some useful formulas that we use through the paper. In Appendix 



we estimate the effect of the correction to the Schwinger-Dyson equation that comes from the non-perturbative 
contribution to the vertex function. At last, in Appendix |C|, we present the approximate analytical solutions to the 
BS equations. 



II. THE MODEL AND NOTATION 



In the case of two flavor dense QCD, the original gauge symmetry SU{i)c breaks down to the SU{2)c by the Higgs 
mechanism. The flavor SU{2)l x SU{2)ii group remains intact. The appropriate order parameter is given by the 
vacuum expectation value of the diquark (antidiquark) field that is an antitriplet (triplet) in color and a singlet in 
flavor. Without loss of generality, we assume that the order parameter points in the third direction of the color space. 



^=(0|£^^£3a.(*D):7M*^)JO>> (2) 



where ^/j and ^£ — C^|^ are the Dirac and its charge conjugate spinors, and C is a unitary matrix that satisfies 
C^^'y^C = — 7J and C = — C"^. Here and in what follows, we explicitly display the flavor {i,j = 1,2) and color 
(a, 5 = 1,2) indices of the spinor fields. It is also appropriate to mention that the subscript and superscript indices 
correspond to complex conjugate representations. 

The order parameter in Eq. (||) is even under parity. Such a choice is dictated by the instanton induced interactions 
[0,|| which, despite being vanishingly small at large chemical potential, could be sufficiently strong for picking up 
the right vacuum. In addition, any bare Dirac masses of quarks (which are non-zero in nature) should also favor the 
parity-even condensate [g^,|^ . 

With the choice of the order parameter orientation as in Eq. (0) , it is very convenient to introduce the following 
Majorana spinors: 



n = V^:+£3ab£*'(V'^)5, a =1,2, (3) 

^L^-^L-esafcE^'C^^)?, a =1,2, (4) 

which are built of the Wcyl spinors of the first two colors, 

^l=V+{-^n)l, (^%^=7'-(*g)5, (5) 

ci^a^V-{-^D)\, (0%^=7'+(vI/g)5. (6) 

Here 'P± — (1 ± 7''')/2 are the left- and right-handed projectors. The new spinors in Eqs. (||) and (|4|), as is easy to 
check from their definition, satisfy the following generalized Majorana conditions: 

i^^^^e'^^'e^.^l (7) 

m: = -e"^e.,K (8) 

In the color superconducting phase of QCD in which quarks are known to acquire a dynamical (Majorana) mass, the 
use of the four-component Majorana spinors, built of the Weyl spinors, is most natural. Of course, when quarks are 
massive and the chiral symmetry is explicitly broken, it would be more appropriate to consider the eight-component 
Majorana spinors, made of Dirac ones. 

With our choice of the order parameter that points in the third direction of the color space, only the quarks of the 
first two colors take part in the condensation. The quarks of the third color do not participate in the color condensate. 
It is more convenient, therefore, to use the left and right Weyl spinors, 

^^=V+{^d)1, {i^^),=V-{^%)l (9) 

cj^^=V-{-^D)l, (0^),=n(^g)?, (10) 

for their description. Notice that the color index "3" is omitted in the definition of V'* and (/)*. 

In the color superconducting phase with a parity even condensate (0), parity is a good symmetry. Then, all the 
quantum states of the Hilbert space, including those in the diquark channel, could be chosen so that they are either 
parity-even or parity-odd. In order to construct such states explicitly, we would need to know the following parity 
transformation properties of the spinors: 

^\x) ^ i°cj,\x'), ^\x) ^ 7V(x'), (11) 

^<^{x) ^ -7Vf (^'), '/'f (^) - -7Vf (^'), (12) 

^'a{^)^l'K{A. $U^)->7"*L(^')- (13) 

where x = {xq,x) and x' — (.tq, ~x). 
Before concluding this section, let us rewrite the order parameter (0) in terms of the Majorana spinors, 

^ = -{Q\^V-K + ^'i'P+Kn. (14) 

This representation is explicitly SU{2)l x SU{2)ji x SU{2)c invariant, and so it is very convenient. By making use 
of the transformation properties in Eq. (|l3|), we also easily check that ip is even under parity, as it should be. 

III. SCHWINGER-DYSON EQUATION 

In order to have a self-contained discussion, in this section we briefly review the Schwinger-Dyson (SD) equation 
using our new notations. This would also serve us as a convenient reference point when we discuss more complicated 
Bethe-Salpeter (BS) equations. 

To start with, let us introduce the following multi-component spinor: 

(15) 

built of the left fields alone. Similarly, we could introduce a multi-spinor made of the right fields. In our analysis, 
restricted only to the (hard dense loop improved) rainbow approximation, the left and right sectors of the theory 




completely decouple. Then, without loss of generality, it is sufficient to study the SD equation only in one of the 
sectors. 

With the notation in Eq. (nS) , the inverse full propagator of quarks takes a particularly simple block-diagonal form, 

g;1 = diag [Sp'Sj'S^^ , s;Mv , s-;^/) , (16) 

where, upon neglecting the wave functions renormalization effects of quarks |lfl-llq], 

Sp' - -* (^ + M7%' + ApT'- + ApV+) 

= -* [{Po S ) 7°A+ + (A;)*A+ + {po + e+) 7° a; + (A+)*A;] V+ 

-* [(Po - <) 7"A+ + A+A+ + (po + e;) 7° A; + A^A^ ] V^, (17) 

Sp' = -^ (:^ + M7°) V+ = -if [{po - e-) A+ + (po + e+) A^] V+, (18) 

Sp^ = -^ {^-t^l°) V- = -if [(po - 4) A+ + (po + e-) K-] V-, (19) 



with e^ — \p\± ^. The notation for the gap function, Ap — A+A+ + A^ A^ and Ap = fAj^f, as well as the 
"on-shell" projectors of quarks 



' ^\ \P\ 



A- = - 1±^ , a = fj, (20) 



are the same as in Ref. fO] . 

Now, it is straightforward to derive the matrix form of the SD equation, 

G;1 = {G^y' +Anas J ^j^'^G,T^''iq,p)V^4q~p), (21) 

where j"^^ and T"^^ are the bare and the full vertices, respectively. This equation is diagrammatically presented in 
Fig. n. The thin and bold solid lines correspond to the bare and full quark propagators, respectively. The wavy line 
stands for the full gluon propagator. 



The only complication of using the multi-component spinor (15) appears due to a more involved structure of the 



quark-gluon interaction vertex. Indeed, the explicit form of the bare vertex reads 

/ [{T^)J> - 26i{T^),'V-] 5^^ {T^yV+6^, -i^^T^yV- \ 
l^^=f\ {T%'V+6^, {T%'V+6^, , (22) 

V -{T^)J^et^V. -{T^)-^V.S,' J 

where £^{, = e'^^e^ac and T^ are the SU{3)c generators in the fundamental representation [tr(r'^r^) — S^-^ /2]. The 
relatively complicated structure of the bare vertex might suggest that our notation are somewhat unnatural. As 
we shall see in Sec. E^, because of the breakdown of the SU{3)c symmetry, this structure, on the contrary, is quite 
natural, and it is especially so in the case of the full quark-gluon vertex function. 

The gluon propagator in the SD equation is the same as in Ref. ||ll|. When the Meissner effect is neglected, the 
propagator in the Euclidean space (fcg = ik^) reads 

V^^^iih, k) ^ S^^V^^ik^, k) c iS^^^ ^ ow 

'^ \k\3 + TrAP\k4/2 ^ 

+,5AB J: O^^J+iS^^ ^O^^J, (23) 

kl + \k\^ + 2M^ ^" kl + \k\^ ^" 

where M^ = N ^asp? ji: (with N j — 2), and Olil are the projection operators of three different types of gluons 
(magnetic, electric and longitudinal, respectively), see Ref. [O. The Meissner effect could be qualitatively taken into 
account by the following replacement of the magnetic term |11[| : 

i^^B W q{\) _^ ^^AB W q{1) /24) 

|^|3-f^Af2|A:4|/2 ^" |fc|3-|-7rM2(|fc4| + c|Ao|)/2 "' 



in the propagators of those five gluons that correspond to the broken color generators {A, B = 4, . . . , 8). In this last 
expression, c = 0(1) is a constant of order one. 

By inverting the expression in Eq. (Jl^), we obtain the following representation for the quark propagator: 



Gn — diat 



:{sX^',^ sp5',, hk'), (25) 






l°{Po - S ) " Ap . 7"(P0 + S)-(A -)* 

>2_(e-)2_|A-|2 P ++%2_(,-)2_|A- 



+ *^2_/,-^2_|A-|2^V' + + * ^2_r.-^2_|A-|2 ^V ' "' ^^'^i 



Po + 4 Po - ep ' 



(27) 



T^A+P-i. "/^A'P^ 
s, ^ Cs!:,Ct = zl^Vj: + ,1^lZ±. (28) 

PO + Ep PO - Cp 



The bare propagator in Eq. ( |2l| ) is similar but with zero value of the gap. 

In the improved ladder (rainbow) approximation, both vertices in the SD equation are bare. By making use of the 
propagators in Eqs. (p3) and (Eq), along with the vertex in Eq. (E2|), we derive the well known gap equation [nTj-nq, 



3'—J (2^)4 52_(^-)2„|A-|2 



Ap = o^as I ,J,, J ;'_,.," \J,i v,Aq-p)- (29) 



After calculating the trace and performing the angular integration (see Appendix |A|), this equation considerably 
simplifies. Then, by assuming that the dependence of the gap on the spatial component of the momentum is irrelevant 
in the vicinity of the Fermi surface, one arrives at the following approximate equation [nOkhsl : 



._, ^ 2as f dq4A (94) A 
A {p4)~-— =ln- r, (30) 

StT Jo Jr,2^\^-\2 \P4-q4 



'91 + 1^0 
where A = {4:TT)'^^^fi/as . The analytical solution to this equation is relatively easy to obtain |ll[ 



A-ip,)^\A^\Joi^uJ^J, P4<|Ao-|, (31a) 

A-(p4) ^ |Ao-|^j2(^) + J2(j.)sin (^In^) , P4 > |Ao-|, (31b) 



where Ji{z) are the Bessel functions and v = y^8a^j9TT. The corresponding result for the value of the gap reads 



,^_, (47r)3/2e/. / 3^3/2 

1^0 I ^ J5/2 ^ ^P ("WV^ ' ' ^^^^ 



where e = 2.718 . . .. Most of the existing studies |10|-|16| seem to agree upon the dependence of this result on the 
coupling constant. The issue of the overall constant factor, however, is still not settled down. The analysis of 
Rcf. psl, for example, suggests that the wave function renormalization effects of quarks give an extra factor of order 
oncR. Another source of corrections might be due to the running of the coupling constant [Bl|. In addition to these, 
as we argue in Appendix ^ using the Ward identities, there also exists at least one non-perturbative correction that 
could modify the constant factor in the expression for the gap. 



^Note, that the argument of Ref. |l3| is somewhat incomplete, since the calculation is performed for the critical temperature 
rather than the order parameter itself. The celebrated BCS relation between the critical temperature and the gap might not 
be satisfied after the Meissner effect is carefully taken into account. 



IV. WARD IDENTITIES 

As in any other gauge theory, in order to preserve the gauge invariance in QCD, one has to make sure that some 
exact relations (Ward identities) between Green functions are satisfied. In this section, we consider the simplest Ward 
identities that relate the vertex functions and the quark propagators. In addition to establishing the longitudinal part 
of the full vertex function, these identities will play a very important role in our analysis of the BS equations for the 
NG and pseudo-NG bosons. 

In general, the structure of Ward identities in non-Abelian gauge theories (the Slavnov- Taylor identities) is much 
more complicated than in Abelian ones: they include contributions of the Faddeev-Popov ghosts. Fortunately in the 
(hard dense loop improved) ladder approximation, used in this paper, the situation simplifies. Indeed, since the direct 
interactions between giuons are neglected in this approximation, the Ward identities have an Abelian-like structure. 

To start with, let us rewrite the conserved currents (related to the color symmetry) in terms of the Majorana 
fields, defined in Eqs. (||) and (^, and the Weyl spinors of the third color. By making use of their definitions, it is 
straightforward to obtain the following representation for the currents: 



+ -l=58^$«(x)7^7'_$;(x) + Ux)^^.T^^'P-<^\{x) + ^U^)V+i^T^^\x) + U^)^^.T^^'t^\^), (33) 

where, in agreement with the remark made above, the irrelevant (in this approximation) contributions of giuons and 
ghosts are omitted. Here and in the rest of this section, we assume that A = 4, ... ,8, i.e., we do not consider the 
currents which correspond to the generators of the unbroken SU{2)c subgroup. 

As we mentioned earlier, we are interested in the Ward identities that relate the quark-gluon vertices to the 
propagators of quarks. Therefore, let us introduce the following (non-amputated) vertex functions: 

^tfj^^y) = (0|rj;f(0)*L(x)V;,(y)|0), (34a) 

Tf:;\x,y) = {0\Tj^iOWix)^^iym, (34b) 

rf';{x,y) = {0\Tj^iO)^lix)^^iym, (34c) 

ttrix,y) = (O|rj;f(O)0^(x)$,"(y)|O), (34e) 

Besides the operator of the conserved current, the first three vertices include only the left-handed quark fields, while 
the other three vertices contain only the right-handed fields. Because of the invariance under the parity, all the mixed, 
left-right vertices are trivial. For this reason, they are of no special interest here, and we do not consider them. 

As usual, in order to derive the Ward identities, one needs to know the transformation properties of the quark fields. 
The color symmetry transformations of the Dirac spinors are well known. By making use of them, it is straightforward 
to derive the following infinitesimal transformations for the spinors of interest: 

(35) 
(36) 

(37) 

(38) 
for the left-handed fields, and 



S^j' 


^zu;^{TtV+K+T,^'^'), 


Si:^ 


= -iLU^ i^fV-T^'' + i^^Tf) , 


SK 


\ ^v3 


s^^ 


\ 2V3 



ILO 



T^''V-<Pl + T^^(j>'), (39) 



— lOJ 



{^^r+T^' + ^,T^^) , (40) 



5^f 



•^ I ^^V+T^" 



l^r^Aa _ ,3ab^^^^Cjj.A3 _ Si^<tfV^ 



(42) 



for the right-handed fields. In aU these expressions, w are small parameters, parameterizing the transformations of 
the SU{3)c group. 

In a standard way, by making use of the current conservation as well as the definition of the vertices in Eq. (|3J) , 
we obtain the following Ward identities for the non-amputated vertices: 



■^A.i 



P^'K^^^^k + P,k)= iT^'5] [sk - Sk+P] V^, 
P^Tf::'{k + P,k)^ iTt5\P+ [Sk - su+p] , 






P'^r^^;;(fc + P,k) = —rft^] [W+Sk - V-Sk - 2Sk+pV- + Sk+pV+] , 



p^K^^ + p,k) = ir^'^ 



;-f^A 



P>'tf-;{k + P,k)^^6i6' 



Sk - Sk + P V+y 

Sk — Sk+P , 
2r-Sk - V+Sk - 2Sk+pV+ + Sk+pV- 



(43a) 
(43b) 

(43c) 

(43d) 
(43e) 
(43f) 



where 5*^, Sk and Sk, Sk are the Fourier transforms of the quark propagators in the left and right sectors, respectively, 



S{x - y)6^d} 
s{x - y)5] 



{Q\TK{x)mvm, 



S{x-y)5'i5]^{Q\T^l{x)mym, 
~s{x-y)6) = {Q\Tcl,\x)Uvm- 



(44) 
(45) 
(46) 
(47) 



As we discussed in Sec. Ill, in the approximation with no wave function renormalization effects, the explicit form of 
the momentum space propagators for the left-handed fields is given in Eqs. (^) and (|2^). For the completeness of 
our presentation, we also mention that the right-handed propagators are the same, except that the projectors V- and 
V+ interchange. 

At this point, let us note that the use of the non-amputated vertices in this section is not accidental. In fact, it 
is crucial for a quick derivation of the Ward identities. Other than that, the non-amputated vertices are not very 
convenient to work with. In fact, it is amputated rather than non-amputated vertices that are usually used in the 
Feynman diagrams. For example, both the bare and full vertices in Fig. n^ are the amputated ones. Similarly, it is the 
amputated vertices that appear in the BS equation in Sec. M. The formal definitions of the amputated vertices read 






-p,k) 
-p,k) 
-p,k) 
-p,k) 
-p,k) 
-p,k) 



Sklp'^aUk + P,k)s 



„-l 



~\A.ai 



CI— 1 j^A.i 
^k+P'- j,ti 

^ ^k+P^aj,tiik 



(k- 

(fc^ 



k ' 

VP,k)S^\ 
P,k)S^\ 

^P,k)rs-\ 



k+p'- ],fi 



s-1 



^k+p^j,'t,(k 



{k + P, k)S 
P,k)S: 



k > 



(48a) 
(48b) 
(48c) 
(48d) 
(48e) 
(48f) 



These, as is clear from our discussion above, are directly related to the quark-gluon interactions. As is clear from 
Eq. (Es), they satisfy the following identities of their own: 



P''Kj:,ik + p,k) = n^f^^ [s^Ip - s-^] n, 

P^vf^^ik + P,k) = ^Tt5YP- [slip S^'] ' 



p^Tf::{k + p,k) 



-^A.i 
aj,fj.\ 



p^'r^f.,{k + p,k) 



2\/3 



Si6^^[2S-lpV+-S-lpV. 



^k+P 



— 1 

Si. 



■P- 



2v-s,'+r+s-'], 



(49a) 
(49b) 

(49c) 
(49d) 



^A.aifj I D ;^\ -rpAa H'', 



P^T^::' {k + P,k)^ iTr5';P+ ~s-' -SZ' , (49e) 



P^^Tf:i^{k + P,k)^^^5t5] 2~S^lpV--S^lpV+-2V+~S^'+V-~S^' . (49f) 



In the rest of the paper, we are going to use these Ward identities a number times. Because of a relatively simple 
structure of the inverse quark propagators, this last form of the identities will be particularly convenient. 

In connection with the Ward identities, it is appropriate to mention here the complementary analysis of Ref. [ p2[ . 
The authors of that paper consider the contribution to the Ward identity that is directly related to the wave function 
renormalization of quarks. 

V. DERIVATION OF THE BS EQUATION 

In quantum field theory, bound states and resonances reveal themselves through the appearance of poles in Green 
functions. These latter satisfy some general BS equations which usually are rather complicated. 

To consider the problem of diquark bound states in cold dense QCD, one has to introduce a four-point Green 
function that describes the two particle scattering in the diquark channel of interest. The residue at the pole of the 
Green function is related to the BS wave function of the bound state. By starting from the (inhomogeneous) BS 
equation for the four-point Green function, it is straightforward to derive the so-called homogeneous BS equation for 
the wave function. 

In the problem at hand, we could construct quite a few different diquark states. Not all of them could actually be 
bound states. For example, one would not expect from a two particle state to form a bound state unless there is some 
attraction in the corresponding channel. Now, in the dense QCD, the dominant interaction between quarks is given 
by the one-gluon exchange. As we know, this interaction is attractive only in the antisymmetric diquark channels. 
Therefore, without loss of generality, it is sufficient to consider only the following bound states: 

{2nf5\p+-p^ - P)x^t\p,P) = (27r)V(p+ -p_ - P)5lx{p,P) 

^{Q\T^l{p+)i,,{-p,)\P-h)L. 6-1,2, (50a) 

{2^)H\p+-p^ P)X^l-,{p,P) - {2nrS\p+-p^ - P)5lKp,P) 

^{Q\Ti;'{p+)^(-p-)\P■ra)L, a =1,2, (50b) 

{2^)H\p+-p^ - P)v{p,P) = {0\T^l{p+m{-p-)\P)L, (50c) 

{27r)^6\p+ -p--P)a{p,P) = {Q\T^'\p+)4,,{-p^)\P)L, (50d) 

plus the states made out of the right handed fields, 

{2nY5\p+ -P-- P)xS\p. P) = {■2tt)H\p+ - p_ - P)5lx{p, P) 

^{Q\T^\{p+)U-P-)\P\h)R. 6=1,2, (51a) 
{2^)HHp+ -p--P)\lj{p,P) = {27TrS\p+-p^-P)Sf\{p,P) 

= {0\Tfip+)^''^{-p^)\P;&)R, 5=1,2, (51b) 

i27rf S\p+ - p_ - P)f,{p, P) = (0|T$l(p+)$?(-p_)|F)fi, (51c) 

{2'Kf5Hp+ -P-- P)^{P, P) = mcj>^{p+)U-P-)\P)R, (51d) 

wherep = (p++p_)/2 and the quantities on the right hand side of these equations are defined as the Fourier transforms 
of the corresponding BS wave functions in the coordinate space. Notice that the analogous states containing the charge 
conjugate fields of the third color, ipf and 0p, are not independent. Because of the property in Eqs. (H) and (H), they 
are related to those already introduced. 

For completeness, let us note that the only other diquark channel that we do not consider here is a triplet under 
SU{2)c. It is however clear that the repulsion dominates in such a channel because this triplet comes from the SU{3)c 
sextet. Notice that, although one does not expect the appearance of a (T bound state, we keep the cr wave function in 
the analysis. This is because the equations for the BS wave functions of the two singlets, or and rj, may not decouple. 
Notice also that the doublet, antidoublet and singlets coming from the S'C/(3)c triplet and antitriplet can mix with 
the doublet, antidoublet and singlets coming from the SU{3)c nonet (octet plus singlet). 

Before proceeding further with the analysis of the bound states, let us recall that parity is not broken in dense QCD 
with two flavors, see Sec. ||. Then, all the bound states can be chosen in such a way that they are either parity-even 



\P;n)p = ^[\P;n)L-\P;n)n], (53) 



or parity-odd. Clearly, the states in Eqs. ( pO| ) and ( pT] ) do not share this property. In order to fix this, we could have 
constructed the following scalars and pseudoscalars, 

\P;n),^^(^\P:n)L + \P;n)j,'^, (52) 

1 

where n denotes the appropriate state. 

In our analysis, however, we find it more convenient to work with the bound states constructed of either left-handed 
or right-handed fields separately. This is because, in the (hard dense loop improved) ladder approximation, the two 
sectors of the theory stay completely decoupled. Besides that, the dynamics of the left and right fields are identical 
in the approximation used. Under these conditions, the degeneracy of the left and right sectors is equivalent to the 
degeneracy of the parity-even and parity-odd ones. In this way, we reveal the parity doubling property of the spectrum 
of bound states in QCD at asymptotically high density of quark matter Jj 

A. Equations for the non-amputated wave functions 

In order to derive the BS equations, we use the method developed in Ref. fsS] for the case of zero chemical potential 
(for a review, see Ref. |34|). To this end, we need to know the quark propagators and the quark-gluon interactions in 



the color superconducting phase. From the analysis of the SD equation (see Sec. Ill), we got the structure of the quark 
propagator. We also know that the approximation with no wave function renormalization effects is quite reliable, 
at least in the leading order. By combining these facts together, wc arrive at the following effective Lagrangian of 
quarks: 



+0, (^ + ^^7") v-^^ + 0^'^ {^ - M7") v+^f + ^. ^s'T^-f - 0^' /s'n^f + *" A'r-<f>' - $^£^, 43V+^f, (54) 

where, by definition, A = A+A+ -I- A~A~, A = j'^A^j^, and e°^ = e^°'^eij. The choice of A, as is easy to check, 
corresponds to the case of the parity-even Majorana mass. 

The effective Lagrangian in Eq. (|5J) is the starting point in derivation of the BS equations for the wave functions 
introduced in Eqs. ( pO| ) and (|5l|). While using the notation of the multicomponcnt spinor in Eq. (llH), it is natural to 
combine the (left-handed) wave functions of the bound states into the following matrix: 



X{p,P) 



' iry(p, P)b^8'^ xi'^ (P, PW, i^cC {\l^ {-p, P)f C^^ 



\^)ip,P)S^, (t(p,P)5\ 

1^ C {x^\~P, P))" C^ef^ Ca^i-P, P)CH. 



(55) 



/ 



where we took into account the property of Majorana spinors given in Eq. (0). We could also introduce a similar matrix 
wave function for the right-handed fields. Since, however, in the (hard dense loop improved) ladder approximation, 
the left-handed and right-handed sectors decouple, we study one of them in detail, and only occasionally refer to the 
other. 

In the (hard dense loop improved) ladder approximation, the BS wave function in Eq. (pq) satisfies the following 
matrix equation: 



G-' [p + I) X(p; P)G-^ (^ " f ) = "^™^ / (04^^''^^«' P^^'^'^'^t^^'i ~ P)' 



(56) 



■^Notice that there are some subtleties in applying this parity doubling argument to the case of the (pseudo-) NG bound states, 
see Sec. 0. 



where V^^ [q — p) is the gluon propagator and 7"^*^ is the bare quark-gluon vertex. This approximation has the 
same status as the rainbow approximation in the SD equation. It assumes that the couphng constant is weak, and 
the leading perturbative expression for the kernel of the BS equation adequately represents the quark interactions. 
Schematically, the BS equation (^^ is shown in Fig. g. 

By writing it in components, we arrive at the set of four equations. 



S;lp/^X^}^{j>,P)s-^lp,^=-l^a, j ^,r[p-x'i\q,P)V-+W-C[x^^^ 



3 'J (2^)4 



''p+P/2 



\l){p,P)s;l 



P/2 



2 

-— ttq:, 
3 






[7^ 



V+X^i){q,P)V+ + 3V+C(\i^{-q,P)) C^V+ 



-V+\Uq,P)V^ YV,,,{q~p) 



(57) 



(58) 



s-lpi2'n{p^P)s^- 



'p+P/2 



p-P/2 



3 



(2^ 



V+in{q, P)V+ + V-ri{q, P)V- + -V-r^iq, P)V+ 



-\v+ri{q,P)V^ + ^V+CT{q,P)V^ + ^r-Ca^i-q,P)C^V+ 



-f'^V^^iq-p), (59) 



s;+P/2'^b. PK-p/2 = - 3™« J (^7^ [3V+v{q, P)'P- + 2V+a{q, P)V-] i^V^^q - P)- 



(60) 



The right-handed fields satisfy a similar set of equations. 

In order to solve the BS equations, it is important to determine the Dirac structure of the BS wave function. There 
is the following useful statement. Let us consider a BS wave function of an arbitrary bound state for a non-zero 
chemical potential in the center of mass frame. Then, the number of independent terms in its decomposition over the 
Dirac matrices coincides with the number of the terms in the decomposition of the BS wave function at zero chemical 
potential. 

The proof of this statement is simple. The Dirac decomposition is determined by all the space-time tensors 
characterizing the bound state, e.g., the momenta P'^ and p'^ , the polarization vector e'^ (in the case of a massive spin 
one bound state), etc. In this respect, the case of a non-zero chemical potential is distinguished by the occurrence 
of only one additional vector u'^ = (1, 0). But in the center of mass frame, where the total momentum P'^ — (Pq, 0), 
the vector u^ is proportional to P"^ and therefore is not independent. Thus, the number of terms in the Dirac 
decomposition of a BS wave function in this frame is the same for both zero and non-zero chemical potential. 

Of course, there is the essential difference between these two cases: while for zero chemical potential, the number 
of the Dirac structures is the same in all the frames, in the case of a non-zero chemical potential, it is different for 
P = (the center of mass frame) and for P 7^ (all other frames). For example, as we will see for spin zero diquarks, 
when ji ^ Q, there are four independent terms in the center of mass frame, and there are eight terms in other frames. 

Strictly speaking, this statement is valid only for massive bound states. However, in the case of spin zero bound 
states it is still valid also for massless states (in particular, for NG bosons): the point is that the limit ill —> is 
smooth for the BS wave functions of spin zero states, and P'^ — > is a very useful limit for studying properties of the 
NG bosons. 

In the next two sections we study these BS equations for the diquark states in detail. In order to approach the 
problem, we first need to determine the Dirac structure of the BS wave functions. As we shall see, the Ward identities, 
derived in Sec. IV, are of great help in dealing with this task. Moreover, in the particular case of the (pseudo-) NG 
bosons, the knowledge of the Ward identities is powerful enough to reveal the complete solution to the BS equations. 
We consider this important case in the next section. 



VI. BS EQUATIONS FOR NG AND PSEUDO-NG BOSONS 

In this section we consider the massless bound diquark states. The latter should include the NG and pseudo-NG 
bosons. Before proceeding to the detailed analysis of the BS equations, it is instructive to describe the qualitative 
physical picture in the problem at hand. 

Let us start from a simple observation. As we stressed many times, the QCD dynamics at large chemical potential 
consists of two essentially decoupled and identical (left-handed and right-handed) sectors. Then, as long as it concerns 
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the diquark paring dynamics, no changes would appear in the model if one enlarges the gauge group of QCD from 
SU{3)c to the approximate SU{3)c,l x SU{3)c,r, assuming that the coupHng constants of both gauge groups are 
identical. In the modified theory, the pattern of the symmetry breaking should be SU{3)c,l x SU{3)c,r — > SU{2)c.l x 
SU{2)c.R- In this case, ten NG bosons should appear. If the gauge group were SU{3)c,l x SU{3)c,r, all the ten 
NG bosons would be unphysical because of the Higgs mechanism. However, since the true gauge group of QCD is 
vector-like SU{3)c, only five NG bosons (scalars) are removed from the spectrum of physical particles by the Higgs 
mechanism. The other five NG bosons (pseudoscalars) should remain in the spectrum. In the complete theory, these 
latter are the pseudo-NG bosons. They should get non-zero masses due to higher order corrections that are beyond 
the improved ladder approximation (an example of such corrections is the box diagram in the BS kernel with two 
intermediate gluons). At the same time, since the theory is weakly coupled at large chemical potential, it is natural 
to expect that the masses of the pseudo-NG bosons are small even compared to the value of the dynamical quark 
mass ]35| |. 

For the completeness of our discussion, let us also add that, even though the massless scalars are removed from the 
physical spectrum, they exist in the theory as some kind of "ghosts" [Q. In fact, one cannot completely get rid of 
them, unless a unitary gauge is found.p It is also important to mention that these ghosts play a very important role 
in getting rid of unphysical poles from on-shell scattering amplitudes [p6[ . 



A. The structure of the BS wave functions of the (pseudo-) NG bosons 



Earlier we mentioned in passing that the use of the Ward identities is crucial for revealing the Dirac structure of 
the BS wave functions of the (pseudo-) NG bosons. Now let us elaborate this point. We start with the definition of 
the vertices in Eq. (|34|). By making use of them, one can show that the corresponding Fourier transforms develop 
poles whenever the total momentum of the incoming quarks, P, satisfies the on-shell condition of a bound state. In 
particular, as P ^ 0, we obtain 



Ttf^,{p + Pl2,p^ 


-P/2) 


Tf;;\p + P/2,p- 


-P/2) 


Tf;;ip + P/2,p- 


-P/2) 


i'tfAp + m,?- 


-P/2) 


ff;;\p + P/2,p^ 


-P/2) 


tf;;ip + P/2,p- 


-P/2) 



P^O 



P^O 



P^O 



P^O 



P^O 



P^O 



^ ^i' 5. 



P^^F^ 

pupix) 



5)T^'x(p,Q), 



(A) 



Y,5Yrtxi,,){p,Q) 



p^'Pi 

7;;]^2^^^«"^(^'0)' 

p(x) p(x) 



(x) 



P^Pi 

pwpw 



w 



J2S}Tf^~Xyp,0) 



P^P^ 
" -i<5;<58^^feO), 



pWp{\) 



pMp(x) 

pup^X} 

-<5;.T3^'^A(p,0), 



P'^P^ 



w 



p-'pi) 



M 2 ' 



(61a) 

(61b) 

(61c) 
(61d) 

(61e) 

(61f) 



where P*-^' and P*^^^ (with x being A, x or v) ^-re the decay constant s of t he (pseudo-) NG bosons. The rigorous 
definition and the calculation of their values will be given in Subsec. VIE. For our purposes here, it is sufficient 
to know that they are constants expressed through the parameters of the theory. Since the Lorentz symmetry is 
explicitly broken by the chemical potential, the dispersion relations of the (pseudo-) NG bosons respect only the 
spatial rotation symmetry. In order to take this into account, we introduced the following four-momentum notation: 
Pff = (P°, —c^P) where c^. < 1 is the velocity of the appropriate (pseudo-) NG boson. 

By recalling that the parity is preserved in the color superconducting phase of dense QCD, we conclude that 
the decay constants of the left- and right-handed composites should be equal, i.e., P^^^ = F^^\ P'^:) = ^(x) ^nd 
pi-n) — p(''). 



^Note that, because of the composite (diquark) nature of the order parameter in color superconducting phase of dense QCD, 
it does not seem to be straightforward to define and to use the unitary gauge there. 
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The existence of poles in the full vertex functions as P ^ is also required by the Ward identities, discussed in 
Sec. [V. Moreover, the Ward identities alone allow us to establish the explicit form of the poles. Indeed, by making 
use of the relations for the amputated vertices in Eq. (Q) as well as the explicit form of the quark propagators in 
Eqs. (H), @ and (H), we obtain 





Ak + P,k) 


T^A.a 
J,fJ- 


\k + P,k) 




\k + P,k) 


f^A.i 


Ak + P:k) 


T^A,a 


\k + P,k) 




\k + P,k) 



p 



ix) 



P^O 



P^O 



P^O 



P^O 



P^O 



_P— 



P^F^ 



(x) 



T^HM^i^V 



ki-'+, 



P 



W 



P'^Pi 



w 



T^''5]l^kV-, 



P. 



iv) 



P'^P'u 



(n) 



S^S' 



V3rx 



AkV+ ~ AkV- 



P 



(x) 



P^P^ 



(x) 



T^'6)AkV-, 



P 



W 



P'^Pb 



m 



T^-S}A„r+, 



p 



in) 



<5«^5! — 



^pin)^»"i 2 



P-'Pb 



AkV- - AkV+ 



(62a) 
(62b) 
(62c) 
(62d) 
(62e) 
(62f) 



Now, by taking into account the definition of the amputated vertices in Eq. (Uq) and comparing the pole residues 
in Eq. (61) with those in Eq. (p^), we unambiguously deduce the Dirac structure of the amputated (as well as 
non-amputated) BS wave functions of the (pseudo-) NG bosons, 



X{P,0) = sz'x{p,0K' ^ 4^v+, 



A(p,0) = SpiA(p,0)(p,0)5p- 



Ap , 



n/3 



rjip,0) = S^'riip,0)Sp' = -^ [Ap 



„V+ - A„V- 



xip,o) = s-'x{p,0K 



Ap 
Fix 

A 



V-, 



A(p,0) = S^'X{p,0)Sp' = -Jw'^+ 



fjip,0) ^ S;'f,{p,0)S;' ^ 



^[ApV^-ApV, 



(63a) 
(63b) 
(63c) 

(63d) 
(63e) 
(63f) 



This concludes our derivation. Before concluding this subsection, we would like to emphasize that the arguments 
used here cannot be generalized for the case of massive diquarks. The reason is that the corresponding on-shell 
pole contributions to the vertex functions [compare with Eq. (pi])] must appear at a non-vanishing momentum P. 
Obviously, the structure of such poles cannot be clarified by utilizing the Ward identities alone. 



B. NG doublet Xa 



(b) 



Now, let us consider the BS equation for the massless x-doublet, see Eq. (p7|). As soon as the color symmetry is 
spontaneously broken in the model at hand, a non-trivial solution to this equation should exist. In order to verify the 
self-consistency of our approach, we have to check that this is the case. 

The most general Dirac structure of the amputated BS wave function, Xa {p, P) = S~}p,2Xa {p, P)^Z-pi2^ ^^^^ ^^ 
allowed by the space-time symmetries is given by 



x'^^\p, P) = Sa^ [xiK + xtK + (Po - S )X^7°A+ + (pn + 4)xtl''K 
+ X3(7 • P) + X4(5 • P) + X5't">"P" + X67V">"P™1 P+, 



(64) 
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where n,m = 1,2,3 are space indices, a""^ = i/2[7",7'"], and the factors (po — Ep ) and {po + e+) are introduced 
here for convenience. It is of great advantage to notice that four out of eight independent functions in this expression 



become irrelevant in the Hmit P — > 0. This agrees with the general statement made in Subsection V A (indeed, there 
are four independent Dirac structures in the BS wave functions of spin zero states at zero chemical potential [p4|). 
We will consider only this limit (which, in the case of NG bosons, implies that the total momentum P ^ 0). 

After multiplying both sides of the BS equation (|57| ) with the appropriate quark propagators on the left and on 
the right, we obtain the equation for the amputated BS wave function. This latter splits into the following set of two 
equations: 

V-x'!^\p,0)V+ = ^Tras J ^j^V+S,x.^^\q,0)s,V-j''V^Aq-p)- (66) 



The quark propagators Sq and Sq are given in Eqs. (26) and (p,T\). As we can see from those explicit representations, 
not all of the terms are equally important. While some of them develop large contributions in the vicinity of the 
Fermi surface, the others are su ppr essed by powers of /i. These latter could be safely neglected in the leading order 
of the theory. Using expression (|2^) for the gluon propagator, we arrive at the following approximate form of the BS 
equations: 

— (xr(-9,0)) ^ _ , ^S^^ 1 A -19 I ^g^-7 T^I^M-P). (67) 




90 + e? ^ ^ 9o - (^9 ) 

V-X^\p.^)V^^~\^aJ ^^^^V^kX ^yV;^ ^_^ X^^{q^Q)i^ll^^^V,M-p)- (68) 

3 ] (27r)* g^ - (e, Y - [A, Y go - eg 

In the component form, these become 

1 C ^4 

xr(p) - X^ois I -p^:;^ [xi{q) + 3xi{-q) - (A~)* (x2"(g) + 3x2"(-g))] 

tr(7'^A^-7"^A;)P^,(g-p), (69) 



3 ,/ (27r) 
1 



■g2_(,-)2_|A-|2 



, -^ -, ^ 1 f dS ['7o-(e^)']x^('Z)-A-xr('Z) 



3 V (27r)4(go_e-)[g2_(,-)2_|A-|2] 
xtr (7'^A|7°7"^A+7°) V,,{q^p), (70) 

plus the expressions that define xt ^^ terms of x7 ■ ^Y noting that the equations for the even and odd combinations 
of the BS wave functions, x7 (p) ^ x7i~P)i decouple and satisfy the same kind of equation, we argue that it is 
sufficient for our purposes to consider only the even combination. Indeed, from the Ward identities, one finds that 
Xi{p) is related to the gap A~. Then, if Xiip) is odd, we must have a non-trivial solution for the gap satisfying 
Aq = Aplp^o = 0. The analysis of SD equation shows that no such solution exists. Therefore, without loss of 
generality, we put x7i~P) ~ X7ip)^ ^^'^ obtain 

-f \ ^ f ^^g XliQ)-{^q)*X2iQ) , I .,_ UK+\T, ( \ ('7^\ 

^^ (^) " 3™^y (^ ql-{e^?-\^-,? ''^'^'' A;)P,.(g-p), (71) 

I -\ -( \ 1 f d^l [go ~(^g')^]x2'(g)-A^xr(g) , f ^,+ ^A+ o^-n f \ f7o^ 

(po-ep)x2(p)--o™« / Tjz^i-f ^vn — r~^i — iX^iTf^'^y^ K^ 7 A+7 )v^,{q-p). (72) 

3 J {2ttY (go - Cp ) 95 - («<? ) - W \ \ 



By comparing our ansatz for the BS wave function in Eq. (p4) with the structure in Eq. ( 33a ) that is required by the 
Ward identities, we see that X2 (p) components should be zero. The direct analysis of the BS equations, on the other 
hand, shows that these component functions X2 cannot be identically zero. It is not hard to pinpoint the origin of 
the discrepancy. Indeed, in our approximation, we completely neglected the wave function renormalization effects of 
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quarks. Upon taking them into account, the Ward identity ( p9| ) would lead to a modified structure of the BS wave 
function, and all allowed Dirac structures would be non-zero. 

Therefore, as in the case of the wave function renormalization, we estimate the effect of X2ip) perturbatively. To 
this end, we use X2ip) = in the leading order of the theory. Then, the equation for Xi{p) reads 

On c omp arison with the gap equation (P9|), we see that Xiip) = i^p)* /^'"^^ [&s required by the Ward identities, see 
Eq. ( |63a| )] is the exact solution to the BS equation in the leading order approximation. Here, of course, we assume 
that Ap is the solution to the gap equation. By substituting the leading order solution Xiip) i'^to Eq. ([72|), we get 
the estimate for X2 (p)- In. the most important region, |Ag"p ^ P4 + (cp )^ ^ M^i "^6 find that 



FWy'pl + («,)= Pl + i^r? 



Now, we also can check that this function could safely be neglected in the equation for Xi (p)- Indeed, its substitution 
into Eq. ( |7l| ) produces the result of order as{A~)* /F^^' which is suppressed by a power of as compared to Xi (p) = 

Therefore, both the corrections due to the wave function renormalization of quarks |10|-^lq] and those due to the 
non-vanishing component functions X2 (p) ^^^ small in the leading order of the theory. Moreover, the consistency 
with Ward identities requires that either both effects are taken into account or neither of them. 

C. NG antidoublet A^g, 

The analysis of the BS equation for the A-antidoublet follows very closely the analysis for the x-doublet. For the 
completeness of the presentation, we still give all the details. 

The most general BS wave function of this antidoublet is given by 

A^,)(p,P) = 5lV- [A+A+ + ArA; + {Po ~ e-)\:,fK+ + (po + e;)A+7°A; 

+ A3(7 • P) + \i{a ■ P) + A5cr">"P™ + A67V">"P™1 . (75) 



As in the case of the x-doublet, to simplify the analysis we restrict ourselves to the case of the vanishing total 

V^,)(p,0)7'-andP_A^,)( 



momentum P ^ 0. Then, the equations for two projections of the wave function, V- A|?~s (p, Q)V- and V- A^- , (p, 0)7^+ , 



read 

74 



P_A^,)(p,0)P_=--7ra,y ^^YV+[s,X\^,-^{q,Q)S, + iS,c(x\~4-q,Q)) C^s^jV+j-'V^Aq - p). (76) 

r^x\^^ip,o)r+ = ^na, J -0^rv+s,xl^{q,o)s,r-Yv^Aq~p)- (77) 

After extracting the most significant in the vicinity of the Fermi surface contributions from the quark propagators Sq 
and Sq, we arrive at 

-D \b I r,\T> 2 f dS ^^ ^+("f°^b ( n^^°(lO-''q)-^q 



3 'J (27r)4' ^ « yq,-eq ('^)^"' ' q^ - (eq)^ - \Aq\^ 

P_A(,,(p,0)P+^--™.y (^T"^^;;^— \a)(<?,0)-^-^-^^-^-^7'-A,7 V,Aq-p). (79) 

Finally, rewriting this in components, we get 
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Kip) 



(Po-ep)A2 (p) ^ -^T^as 



d^q X^{q)-A-X-{q) 



{2nrq',-{e~y-\A-\^ 

A [ql-ie,r]X,{q)~iA-rX^iq) 



{2nY (5o_e-)[52_(,-)2_|A-|2] 



^tr(7^A+7'^A;)l?^,(g-p), (80) 

tr(7^A+7VA+7°)2?p.(?-p), (81) 



where again, without loss of generahty, we assumed that A~ (p) are even functions of momenta. 

By repeating the arguments of the previous subsection, we would find that X2 [p) should be zero in a consistent 
approximation when the wave function renormalizations of quarks are neglected. As in the case of the x-doublet, the 
equa tion for X^{p) component has the solution Xi{p) — —Ap/F'-'^^ that is consistent with the Ward identities, see 
Eq. (HI). 



D. NG singlets r/ 

The case of the massless singlets is very special. This is already seen from the fact that the BS equations for the 
r] and cr-singlets are coupled in general. This might appear somewhat puzzling if one traces back the origin of the 
singlets. While the ry-singlet contains the antisymmetric tensor product of two fundamental representations of SU{S), 
the (T-singlet comes from the product of the fundamental and the anti- fundamental representations. Based on this 
observation, one might have concluded that the bound state should form only in the 77-channel. As we shall see below, 
this argument is not completely groundless, although, the real situation is slightly different. We would like to point 
out that there is no symmetry in the color superconducting phase of dense QCD which could prevent the coupling 
between the two singlet channels. 

Let us consider the equation for the amputated BS wave functions of singlets, ri{p, 0) and a{p, 0) (we consider only 
the case of P ^ below) . The most general structure of the wave functions is 



r,ip, 0) = [77^ A+ + 77+ Ap + [po - Ep )% 7°A+ + (po + ^p )%+7"Ap ] r+ 
a{p, 0) = V-j° [(po - e;)^"A+ + (po + e+)<j+Ap] V+. 



(82) 
(83) 



Before proceeding any further, wc would like to remind the definition of the ry-singlet wave function in Eq. (50c). 
Unlike other diquarks, it is built of only the Majorana spinors. By making use of the generalized Majorana property 
in Eq. (||), we observe that the BS wave function ri{p, P) should satisfy the following constraint:^ 



Cv^{-p,P)C^^rj{p,P). 



(84) 



While rewritten in components, this restriction is satisfied when the odd components iji ^{p) are even functions of 

momenta, and when 7/2 (—p) — vtip)- 

The equations for different chiral projections of the amputated BS wave functions, 'r]{p,0) = Sp^ri(p,0)Sp^ and 
^■(PjO) = s~V(p, 0)s~-^, read 



V+v{p,0)V+ = 


~-na 


V-vip, 0)V- = 


8 
--na 


P_77(p,0)P+ = 


2 

— -na 




r+vip, o)V- = 


2 
~-na 


V-a{p,0)V+ = 


2 
--na 



(iTT)-* 

d-^q 

(2^ 

d^q 
(27r 

d^q 
(27r) 

^J^V+ {3Sgfj{q, 0)Sg + 2s,a{q, Q)s,) V-^^V^^iq - p). 



rV+S,ri{q, 0)S,V+YV^,{q - p), 
^rV+ {^S,ij{q, 0)5, + 6s,a(<7, 0)s,) V-^T^^^M - P), 
'^rV- [5S,r,{q, 0)Sg + 6s, (a(-g, 0))^ s,) V+^-D^^iq - p). 



(85) 
(86) 
(87) 



(89) 



^Strictly speaking, from Eq. (50c) one derives a relation for the non-amputated BS wave function. Assuming that the gap is 
an even function of the momentum, it is straightforward to show that the same relation holds for the amputated wave function 
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In components, these become 



4 f d'g [<;S-(e-)>,-(5) + (a,-)V(9)-2(A,-)-[<;„»-(e-)>,-(5) 



xtr(7^A-7''A+)l?^,(q-p), (90) 

_, , 4 f dU H - (e-)^]r;3-(g) + (A-)^7?r(g) - 2A^'[gg - {e-f]v^{q) 

xtr(7^A+7''A-)l?^,(q-p), (91) 

1__ /■ A (r_ [go - K')^]^2'(g) + |A-pr/4-(g) - A-r;r(g) - (A-)*773-(g) 



(Po ~ Cp )?72 b) - o^"« / TTZh 5 



3 J {^^Y\ [g„2-(e,-)2-|A-p] 

e-)tr (7^A+7"7''A+70) V^M-p). 
H - (£,-)']% (9) + \^-,?Vl{q) - A-77r(9) - (A-)*ry3-(q) 



2^^T^) (90 - ^,")tr (7^A+7"7''A+7°) V,^{q-p), (92) 

J4^ / [„2_ ^^-^l2l„-^„^,4_ |A-|2 

(* - f,-)tr (7''A+7»7''A;7'') !>„„(, - p), (93) 

along with the irrelevant for our analysis expressions for the plus-components. Note that we did not write down the 
equation for rj^{p) component since it is related to 7^2 (p) ^^ ^^ argued above. 

Now, let us analyze the BS equations for the singlets. By repeating the argument involving the Ward identities, 
we see that the component functions ?72 4(-P) should be exactly zero. As oppos ed to the case of (anti-) doublet, a 




crucial difference appears in the case of singlets. As we discussed in Subsec. VIA, the Dirac structure of the BS wave 



function of the ?7-singlet is determined by th e po le structure of the vertex in Eq. (31c) . Th e explicit form of the latter 



is determined by the Ward identity in Eq. (49c), and the result is presented in Eq. (62c). The remarkable property 



of this result is that it does not get any corrections even after the wave function ren orma lization effects of quarks are 



taken into account. To see this, one should note that the mentioned Ward identity (49c) is given in terms of a single 
propagator, Sp. Because of this, all wave function renormalization effects always cancel from leading pole contribution 
to the vertex of interest. 

It is very rewarding, therefore, to check that r]^{p) = —VS{Ap)*/F'^^\ Vsip) = V^Ap/F'-^^ and ri2 4^{p) ~ cr^(p) == 
is the exact solution to the BS equations,H assuming that Ap is the solution to the gap equation. It is interesting 
to notice that no admixture of the cr-singlet appears in this solution for the (pseudo-) NG boson. 

E. Decay constants 

In order to define the decay constants of xi , A(.~>, 77 (pseudo-) NG bosons (as well as their counterparts built of 
the right-handed fields), it is convenient to introduce the following combinations of the currents: 

7 

Jia),^T.^a'j^. «-l,2, (94) 



A=4 
7 



J 



(b) _ 



A=4 
„-8 



J2Tfj^ 6=1,2, (95) 



Jp-J-', (96) 

where j^ for A = 4, . . . , 8 is defined in Eq. (^3|). It is easy to check that the doublet xi couples only to j(a)p, while 
the antidoublet A.^s couples only to j), . The singlet couples only to j'^. 



^We believe that this is the only non-trivial solution to the BS equation, although we were unable to rigorously prove that no 
other solutions exist. 
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We will consider the low energy limit when the energy of the diquark (pseudo-) NG bosons Pq -^ 0. Then, we 
define their decay constants as follows: 

(0|j(a);^(0)|P,6)L=z<5^F(/)F(x), (97) 

{0\jl^\0)\P,&)L^^SiP^,^^F^^\ (98) 

{0\j^mP)L=tPl:^^F^^\ (99) 

where Pff' — (Pqj—c^P] and c^ is the velocity of (pseudo-) NG bosons. 

From the definition of the current in Eq. (p3) in terms of quark fields and from the definition of the BS wave 
functions in Eq. (pOl), we obtain 

pMpix) ^ 1 1 ^tr [7,nx(g, P)]-y (04 tr [j.V+S.+p/Mi, P)s,-p/2] , (100) 

p!.'^F^'^ ^\j (04 tr [V-^,\{q, ^)] = ^ / (0*^ [^-7m^,+p/2 A(g, P)S,^P/2\ , (101) 

P^-^f (") = ^J (0jtr [7,n»7(<Z, P)]-^l (04 tr b.P+S.^PMil, P)S,-P/2] , (102) 

where the doublet and the antidoublet BS wave functions are defined so that Xa = S'^x a-nd A,^s — 6^''X. The generic 
definition is diagrammatically presented in Fig. ^. 

The definitions of the decay constants above are exact. The problem is, however, that the solution for the BS wave 
functions x(g,P), X{q,P) and ri{q,P) at P ^ is very hard to obtain. In order to get the estimates of the decay 
constants and velocities, we will use the analogue of the Pagels-Stokar approximation [B^ (for a review see Ref. p8[). 
In this approximation, the wave functions at P 7^ are substituted by their values at P = 0, i.e., 

x(g,P)~x('Z,0)- A_7'+, (103) 

A(<z,P)c.A(g,0) = --^P_, (104) 

77(g, P) c. f^{q, 0) = -^ (a,P+ - AgV^) , (105) 

where the right hand sides are fixed by the Ward identities. By making use of this approximation and the explicit 
form of the quark propagators in Eqs. (Eq) and (p7|), we derive in the weak coupling limit. 



clPJ IGn^J, J, \PcosHjJ [^| + (,-)2+^|A-|2]2 8n^\^P 

V J\ciPJ 16^3 jg J^ {Pcos^S U [ql + {e^)^+x\A^\^f S^^' U^ J 



m' r.n.:^n( ^0 1 / dq,de-\A-\^ _ fi^ ( P, 



i^'-'n:>}~-6i:-Hpi^.}i 



[g| + (6,-)2 + |A,-p; 



2 8^2 ) Ip 



3 



(108) 



We observe that the decay constants of all NG bosons are of order /i, and all the velocities are equal to l/\/3. This 
agrees with similar results for the NG bosons in three flavor dense QCD pSJ-pq], as well as with the studies of Ref. [ p9[ , 
dealing with the two flavor QCD in the framework of the effective theory approach. Notice that these estimates of 
the decay constants imply that the infrared masses of the five gluons are of order g^/i. (For a discussion of subtleties 
concerning the generation of the mass of the eighth gluon, see Ref. |40||). 

For completeness of our presentation, let us mention that the above expressions for the decay constants in the 
Pagels-Stokar approximation also contain the following subleading derivative term corrections: 



F(x)'^' 



pW^ 



^P f\,, ^^^^ f dq^de\A-\^q^d,,\A- 



, „ , x(l — x)dx I o 

4^'-/n ./ [ql + {e-)^+x\A^r' 
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M 



(f(x)y 



= s 



pW 






dx 



dqideqi{i/4.-x)dq^l^-\' 



127r3 






(109) 




2 /•! 



dq4&|A-|2eae|A-|2 



[g| + (6,-)2 + a:|A,-p]^ 






247r3 



A* 



da; 







dqid€€{2,/A~x)d^\/\-[' 



167r3 7 [5| + (e-)2 + |A-|2]^ 



487r3 7 [g| + (e-)2 + |A,-p] 



2 ' 



(110) 

(111) 

(112) 



where e = e^ = |q| — /i. In calculation, we also assumed that the gap is an even function of 34 and e. By making 
use of the explicit solution to the gap equation (see Sec. Ill and Appendix B in Ref. |11|), we could check that thes e 
corre ctions are suppressed by a power of the coupling constant as compared to the leading order results in Eqs. ( 106 ), 
( |T07| ) and ( |lO^ ). 

Here we worked out the decay constants of the left-handed diquark composites. It is straightforward to repeat the 
analogous analysis for the right-handed diquarks. As it should be, the results would be the same, in agreement with 
the invariance under parity. 

It is worthwhile to note that the decay constants in the two sectors of the theory could have been equivalently defined 
through the matrix elements of the left- and right-handed color current. As we discussed earlier, both such currents 
are approximately conserved at large chemical potential. In other words, the latter means that the corresponding 
vector and the axial vector currents are also approximately conserved. Now, by making use of the definition of the 
scalar and pseudoscalar diquarks in Eqs. (|5^) and (^, it is easy to show that their decay constants are equal, up 
to a factor of v^, to those of the left- and right-handed states. Of course, here we assume that the definition of the 
decay constants of scalars and pseudoscalars are given in terms of the vector and axial vector currents, respectively. 

While the NG scalars (a doublet, an antidoublet and a singlet) are not the physical particles because of the 
Higgs mechanism, the five pseudoscalars remain in the physical spectrum. Since the latter are (nearly) massless 
in dense QCD, they should be relevant degrees of freedom in the infrared dynamics. The decay constants of these 
pseudoscalars are the physical observables that could be measured in an experiment. The most likely decay products 
of these pseudoscalars should be gluons of the unbroken SU{2)c and the massless quarks of the third color (which 
might eventually get a small mass too if another (non-scalar) condensate is generated |7||4l[]). 



VII. BS EQUATION FOR MASSIVE COMPOSITES 

The essential property of the quark pairing dynamics in Nf — 2 dense QCD is the long range interaction mediated 
by the gluons of the magnetic type [p|JTc|]. Of course, the Meissner effect in the color superconducting phase produces 
masses for five out of the total eight magnetic modes. Nevertheless, there are still three modes that remain long 
ranged. This simple feature has many interesting consequences. One of them was conjectured in Ref. [ p7| where it 
was suggested that there should exist an infinite tower of massive radial excitations in the diquark channels with the 
quantum numbers of the (pseudo-) NG bosons. This conclusion was reached by making use of an indirect argument 
based on special properties of the effective potential in the color superconductor. 

In this section, we study the problem of massive radial excitations by using the rigorous approach of the BS equation, 
derived in Sec. |^. As we shall see, the conjecture of Ref. [^ is essentially correct. At the same time, it will turn out 
that some details of the pairing dynamics are rather sensitive to the specific details (such as the Meissner effect) and 
could not have been anticipated, based on the qualitative arguments of Ref. |2^. 



A. Bound states and the Meissner effect 

In this subsection, we would like to clarify the role of the Meissner effect in the dynamics of diquark bound states. 
The main point we want to emphasize is the existence of two classes of bound states, for which the role of the Meissner 
effect is very different. The first class consists of light bound states with the masses M <C |A(^|. The binding energy 
of these states is large (tightly bound states). The second class includes quasiclassical states with the masses close to 
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their threshold fc|Ag | where k = 1 for diquark doublets, built of one massive and one massless quark, and fc = 2 for 
diquark singlets built of two massive quarks. The binding energy of the quasiclassical states is small, i.e.. 

It is clear that the quasiclassical bound states are sensitive to the details of the infrared dynamics. Indeed, for long 
range potentials, the quasiclassical part of the spectrum is almost completely determined by the behavior of the 
potential at large distances. In the particular case of cold dense QCD, the interaction between quarks is long ranged 
in the (imaginary) time direction and essentially short ranged in the spatial ones ||l^,^. Because of that, the region 
with |A:o| < \^o\ ~ 1^1 ^^ particularly important for the pairing dynamics of the quasiclassical diquark states. This 
implies that the inclusion of the Meissner effect is crucial for extracting the properties of the states from this second 
class (see Appendix y). 

On the other hand, the Meissner effect is essentially irrelevant for the light bound states. This point can be 



illustrated by the BS equations for the lightest diquarks, the massless NG bosons. As was shown in Sec. VI, the BS 



equations for them are essentially equivalent to the gap equation. And we know from the experience of solving the 
gap equation [[ll| - [lq| that the most important region of momenta in the equation is given by |Aq | ^ |fco| ^ |A:| < /i. 
In this particular region, the two kernels, with and without the Meissner effect taken into account [compare Eqs. ( p^ ) 
and (|4|)], are practically indistinguishable. Obviously, the same should be true for all very light bound states with 
A/«|Ao|. 

In the rest of this section, we show that, because of the Meissner effect, an infinite tower of (quasiclassical) massive 
diquarks occurs only in the singlet channel. In the doublet channel, the only bound states are the (tightly bound) 
NG bosons. This is connected with the fact that the interaction in the doublet channel is mediated exclusively by the 
gluons subject to the Meissner effect. In the singlet channel, on the other hand, the interaction is partly due to the 
unscreened gluons of the unbroken SU(2)c subgroup and, therefore, the formation of massive (quasiclassical) states 
is not prohibited. 

B. Massive doublet Xa 

Let us start from the analysis of the BS equation for the massive x-doublets. By choosing the spatial component 
of the center of mass momentum of the bound state zero, P — (M^, 0,0,0), we find that, similarly to Eq. (pi]), the 
most general structure of the BS wave function in the center of mass frame reads 

xi^) {p, M^) - 5} [x{K+ + xtK + {Po - S + Mj2)x2l°A^ + {po + < + M^/2)xtl°A;] V+, (114) 

where, for convenience, we introduced here the factors (po ^ £p + -^^x/^) ^^"^ (Po + ^p + -^^x/^)- Iii contrast to the 
case of the (pseudo-) NG bosons, the structure of the wave functions of massive states cannot be established from the 
Ward identities. Although the vertex functions would also have the poles, corresponding to the massive intermediate 
states, the Ward identities are insufficient for extracting the structure of the residues unambiguously. 



Be repeating the analysis similar to that in Subsec. VI B, we arrive at the following set of equations for the 
components of the BS wave function: 



Xi (p) 



4 f A iqo-e-+Mj2)[x^{q)-{A-rx2{q)] 



rTTO:., 



3 ^ J (2^)4 (go _ e- » M^/2) [(go + M^/2)2 - (e-)2 _ |A,-|2] 

xtr [7^A~7''A+] V^^q - p), (115) 

3 J (27r)'' ((Jo -eq - M^/2) [(<?o + M^/2.y - (e, Y - jA, \^\ 
xtr [7'^70A-7''A;70] V^,{q - p). (116) 

The other two equations express Xi2{p) i^ terms of Xi2{p) and, therefore, they are irrelevant. As is easy to check 
by tracing back the derivation of the BS equations, the interaction in this doublet channel is provided exclusively by 
the gluons that are subject to the Meissner effect in the color superconducting phase. This fact could be taken into 
account qualitatively by replacing the propagators of the magnetic modes according to the qualitative rule in Eq. (p3) . 



In accordance with the discussion in Subsec. VII A, while such modification was irrelevant in the gap equation and 
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in the BS equations for the (pseudo-) NG bosons, it is going to play a very important role in the analysis of the 
quasiclassical massive states. 

In order to proceed with the analysis of the BS equation, we will use the approximation with X2 ~ 0- Recall that 
such an approximation was completely justified in the case of the (pseudo-) NG bosons. It is certain that it gives 
a very good approximation for light bound states, with M <^ |Aq |, in general. Our conjecture is that the ansatz 
with X2 ~ yields a reasonable approximation even in the case of quasiclassical states. In order to justify this 
approximation one would need to prove that the (perturbative) correction due to non-zero X2 i^ small. By repeating 
the arguments used for the (pseudo-) NG bosons, we could show again that X2 i^ irrelevant in the ultraviolet region 
I ^0^1 ^Po ^ fJ- While this observation is promising, it is not sufficient yet because the infrared region < po < \^o\ 
is also important for the pairing dynamics of the quasiclassical radial excitations. Therefore, for such states, it should 
be considered as a (reasonable) conjecture. 

Now, we drop X2 from the analysis and assume that the component functions x^depend only on the time component 
of the momentum (compare with the analysis of the gap equation in Refs. pX|-p^). Then, we arrive at the following 
equation: 

2a, f^ . . .,A, . _. .. A 



Xi{P4)^^ dq,K^^Hq4)Xi{q4)\n- . ,._, , (117) 

yTT Jo m -Pi\ + c|Ao I 



where A = (47r)'^/^/i/Q;s and the kernel reads 



M^U4 + V9l + |AoP -|Ao-|' 
^^^^^(94) = , ^ ^ -r. (118) 



g| + |A„-p 4M2g|^(|A„-p-M2) 



To analyze the BS equation, we could convert it into a differential equation, using the same approach as in the case 
of the gap equation in Ref. [|l^ (see Appendix y). It is straightforward to show then that, in the ultraviolet region 
|A^| < p4 < A, the BS wave function of the massive doublet, XiiPi)^ has the same behavior as the BS wave function 
of the NG boson in the x-channel (which, as we know from the Ward identities, is proportional to the gap function). 
The deviations might appear only in the infrared region, < p4 < |Ag"|. Note, however, that the BS wave function 
in that regi on is essentially constant. To see this, one should notice that the dependence on 734 on the right hand side 
of Eq. dllTj ) comes through the argument of the logarithm, in which either large 94 or c|Aq | dominates. 

By matching the infrared and the ultraviolet asymptotes of the solution, we find that the non-trivial wave function 
Xi{P4:) cannot ha ve z eros in the whole range of momenta, < p4 < A. This means that the only solution to the 
integral equation ( |117| ) corresponds to the NG boson with M^ = 0. Indeed, a solution that describes a massive radial 
excitation of the NG boson must have at least one zero. Because of the Meissner effect, no such solutions exist. 

From physical point of view, the massive gluons cannot provide a sufficiently strong interaction to form massive 
radial excitations of the NG bosons in the doublet channel. To understand this better, it might be instructive to 
point out that, in absence of the Meissner effect, an infinite tower of (quasiclassical) massive states would appear 
(see Appendix O). However, the binding energy of all these states would be very small compared to the value of the 
superconducting gap. This indicates that it is the long range dynamics that is primarily responsible for the formation 



of such quasiclassical bound states. Then, in agreement with the discussion in Subsec. VII A, even relatively small 



screening effects for gluons in the infrared region are sufficient to prevent binding of quarks in the doublet channel. 

C. Massive antidoublet A/gj 

The analysis for the A, ~x -antidoublet resembles a lot the analysis in the previous subsection, so we outline only the 
main points, omitting the unnecessary details. 

The general structure of the BS wave function is given 

\l)ip, Ma) - SJ'V- [A+A+ + ArA; + (po - e; + M^/2)\^^'>A+ + (po + e+ + Ah/2)X+^^Ap] . (119) 

The equations for the components of the BS wave functions of the antidoublet are almost the same as those for 
the doublet. The only difference is that M^^ is replaced by —M\. By repeating all the arguments of the previous 
subsection, we again conclude that, after the Meissner effect is taken into consideration, the gluon interaction is not 
strong enough to provide binding of the massive radial excitations in the A^^n antidoublet channel. 
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D. Massive singlet r/ 

Let us consider massive singlet diquark with ikf^ 7^ 0. Since the equation for rj and a do not decouple, the massive 
radial excitation of the 77 NG boson would have a non-zero admixture of cr. The general structures of the BS wave 
functions in the center of mass frame, P = (M^, 0, 0, 0), are 

rjip, Mr,) - [11^ k+ + r/+A; + (po - ^ + Mr,/2)i^^fK+ + {po + e+ + M,,/2)4j'>A;] V+ 

+ [4 a; + %A; + {po + (-- M^/2)r74-7°A; + (po - < - Mr,/2)4j''A;] V-, (120) 

a{p, M^) - 7'_70 [(po -e-+ M^I2)g- A^ + (po + < + M^/2)a+A;] 7'+. (121) 

The components satisfy the following set of equations: 



1 



-..4 (^ 



9o - (e, 
(A-) ^ 



M, 



1) \2 



771 (9) + (A-*)'% (9) 



2 



P0 + — -Sj%(p)-3'^".y (^ 



'Zo^ - (^.- - ^)^ 



d^g 



(%"(9)+'?r(9)) |tr[7^A-7"^A+]l?^,(9-p), (122) 

1 



)4 [(go - M,/2)2 - (,-)2 _ |A,-p] [(50 + M,/2)2 - (,-)2 _ |A,-p] 



M„ 



X ^ 1 90 + ^ + Eg 



90 - K" - ^f ) %"(?) - (A-)*% (g) 



90 



^+6-) (|A-|V(9)- A-ryr(<z)) jtr [7VA+7VA+] P^.(g - p) 



27ras 



% (P) = g^ra. 



2 



(27r)4 go - Af^/2 - e. 



^tr[7VA+7VA+]P^,(g-p), 
1 



(123) 



(27r)4 [(go - M,/2)2 - (e-)2 _ |A,-|2] [(go + M,/2)2 _ [e-f - \^-\^\ 



<& (^,- + '^f 



Mr, 

2 



%('?) + (A-)'?7i(9)-A, 



2 



(90 + ^)^ -(6-)^ 



% (q) 



^a; 



I Po - — + Ep I ??4 (P) == gTra. 






774-(g)Ur[7^A+7''A;]l?^,(g-p), 
1 



(124) 



(27r)4 [(go - M,/2)2 - (e-r - \A^\'] [{qo + M,J2Y - {e^r - \^^\'] 



M„ 



90 



M 
90 - K" - -yf ) V4{q) (A-)*% (9) 



90 



2 ^« 



(|A-|\2-(9) - A-ryr(g)) tr [7VA-7VA;] V.^Q-p) 



2iTas 



M, 



d'^q g (-9) 

(27r)4 go + M^/2 + 

d^g 



-tr[7VA-7VA;]l?^.(g-p), 



(125) 



1 



^"^ 2' '^j'' ^^^ ^"^/ (27r)4[(go-M„/2)2-(,-)2_|A,-p][(go + AV2)2-(,-)2_|A,-pl 



M 



90 



>; 



M 
ql K" - ^)' ) % (9) - (A-)*% (g) 
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90 - ^ + e-) (|A- \^i^l{q) A-,yr(<z)) jtr [7°7^A+7°7''A;] V,,{q - p) 
^nas I S^ ^^ ^tr [7°7^A+707^A;] V^.{q ~ p), (126) 



3 "7 (27r)4 go _ A/^/2 

In the case of massless NG bosons, the component functions 772^4 equal zero. Similarly to the case of doublets, we 
assume that the ansatz with these functions being equal to zero yields a good approximation also for massive diquarks. 
By substituting 773^4 = into the BS equations above, we obtain the following simple set of equations: 



Vi [P) 3^. J (2^)4 [(g^ _ Mj2y - (e,-)2 - |A,-|2] [{q, + M,/2)2 - (e"): 

xtr[7^A-7''A;]l?^,(g-p), (127) 



|A,-P] 



Va (P) = o 



4 f A [go'-(6-+AV2)^J,?3-(g) + (A-)V(9) 



TTttc 



3 ' J (2^)4 [(q, - M,/2)2 - (e-)2 _ |A,-p] [(go + M,/2)2 - (e-)2 _ |A-|2] 
xtr[7^A+7''Ap]l?^,((7~p), (128) 

plus the equation for a~ which does not allow a non-trivial solution for a bound state. 

Now, in order to solve the set of equations for rj^ and 77;^, we make the following substitution: 

Viip)^~^-^h^ip), (129) 

Vsip) = ^'»3(P). (130) 

and, as in the case of doublets, we assume that the wave functions (ft.1,3) depend only on the time component of the 
momentum p4 = ipo (also compare with the analysis of the gap equation in Refs. p0|-p6[). At the end, we arrive at 
the following equation for the BS wave function {hi = h-^) of the massive singlet: 

hl{p4) = -^ dqiK^'^\qi)hi{qi)\n- :, (131) 

47rJo m-Pi\ 

where A = (47r)'^/^ /i/ofs , and the kernel reads 

,» + |A„-p-(M,/2f 

At this point it is appropriate to emphasize that, as we saw already in the previous two subsections, the Meissner 
effect plays an important role in the analysis of the massive bound states. Indeed, our analysis indicates that only 
the long range interaction mediated by the unscreened gl uons of the unbroken SU(2)c is strong enough to produce 



massive bound states. This is taken into account in Eq. (131) where the effective coupling constant differs by the 



factor 9/8 from the coupling in the gap equation (30) 



In order to get the solution for the BS w ave function hi{p), we use the same method as in the case of the gap 



equation |11|. In particular, we convert Eq. (131) into the differential equation 



Pih'lipi) + h[{pi) + ^K('')(p4)/ii(p4) = 0, (133) 

47r 

along with the boundary conditions, 

/i'i(0)=0 and /ii(A) = 0. (134) 



Now, we solve the differential equation ( |133D in the following three qualitatively different regions: < P4 < 

yJAo"J2T"(M^^/2f , y'lAg |2- {Mn/2f < pi < \Aq\ and \A^\ < pa < K. The kernel ( pl ) has a simple behav- 
ior in each region, and the BS equation allows the analytical solutions. 
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hi{p4) = Co Jo 



"slAo b4 



■npi 



^7r[|A-|2-(AV2y 



for 0<P4<\/|A(7P 



(Af,/2) 



"^'^01 I, for .\/|AoP-(A^V2)'<P4<|Ao|, 



/ii(P4)=C3sinfiW— InA) , for |Ao|<p4<A, 

V2 V TT P4/ 



(135a) 



(135b) 



(135c) 



where J„ and iV„ are the Bessel functions of the first and second type. The solutions are chosen so that the boundary 
conditions are automatically satisfied. In the above expressions, Ci (i = 0, . . . , 3) are the integration constants. To 
obtain the spectrum of the massive diquark states, we match the logarithmic derivatives of the appropriate pairs of 

the solutions at -W |Aq p — (A/,,/2) and |A,^|. After taking into account the equation that determines the value of 
the gap (see Appendix B in Ref. Q), 



ln^.^arctan^'^°('^) 



I An 



JM 



Sets 

977 ' 



the matching condition reads 



Jo (2:0)^1(20) 



Ji{zo)No{zo) + Jo{zo)Ni{zo) 
where the coupling is assumed to be small, a^ ^ 1, and 



sTT / 37r \ 



zo 



as I An 



^7ry'|A„P-(AV2)2 



(136) 



(137) 



(138) 



It is straightforward to check that the left hand side of Eq. ( |137D is an oscillating function having an infinite number 
of zeros (zo ~ 2.40, 3.83, 5.52, . . .). In the weakly coupled theory, each zero (or rather a nearby point that approaches 
the zero as a,, — > 0) determines a corresponding value of the diquark mass^. In the quasiclassical limit, i.e., when 



Mr, — > 2|Aq I from below, the left hand side of Eq. (137) is approximately given by cot(2zo)- Then, we derive the 



following simple estimates for the masses of the 7y-singlets: 



Mi 



4|Ao-p 



8„,2 



2«a- 



TT^{2n + lY 



n> 1. 



(139) 



This agrees with the expression presented in Eq. (|l|) whe n k = 1? jn^ « 0.27. Accide ntal ly, one could also check from 
the position of the zeros on the left hand side of Eq. ( |137[ ) that the expression in Eq. (139) gives a good approximation 
even for the low lying states (n = 1, 2, . . .). Notice that the state with n = does not appear. 



VIII. CONCLUSION 

In this paper we studied the problem of diquark bound states in the color superconducting phase of iVj = 2 cold 
dense QCD. We used the conventional method of BS equations that suits the problem best. We derived the general 
BS equations, and then analyzed them in spin zero channels. 

Our analytical analysis of the BS equations in cold dense QCD shows that the theory contains five (nearly) massless 
pseudoscalars (pseudo-NG bosons) which transform as a doublet, an antidoublet and a singlet under the unbroken 
SU{2)c- To the best of our knowledge, these pseudoscalar diquarks have not been discussed in the literature before. 
We estimate the decay constants of these pseudoscalars, and find that their orders of magnitude are the same as that 



^Note that there is also a zero at zo = 0, but Eq. (137) does not have a solution in its vicinity. 
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of the chemical potential. The velocities of the pseudoscalars are equal to l/-\/3, and this coincides with the velocity 
of the NG bosons in three flavor QCD. While being (nearly) massless, the five pseudoscalar diquarks should be the 
relevant degrees of freedom in the low energy dynamics of Nf = 2 dense QCD. 

The parity-even partners of the pseudoscalar diquarks are the NG bosons which are the ghosts in the theory. 
Although they are removed from the spectrum of physical particles by the Higgs mechanism, one cannot get rid of 
them completely, unless a special (unitary) gauge is defined. Since the order parameter is given by a diquark composite 
field, it does not seem to be straightforward to define and to use the unitary gauge in dense QCD. In all the covariant 
gauges we use here, the NG bosons are always present and they play an important role in removing unphysical poles 
from physical scattering amplitudes. 

We also studied the problem of massive diquarks. In accordance with the conjecture of Ref. |2^], there exists an 
infinite tower of massive bound states which are the radial excitations of the (pseudo-) NG bosons. As a result of 
the Meissner effect, it appears that the massive radial excitations occur only in the singlet channel. This could be 
understood in the following way. The interaction in the doublet and the antidoublet channels is provided exclusively by 
the gluons affected by the Meissner effect. Such interaction is not sufficiently strong to form massive radial excitations 
in those channels. The important point in this analysis is the different role the Meissner effect plays for tightly bound 
states and quasiclassical bound states. 

As we know, the parity is unbroken in the color superconducting phase of two flavor dense QCD. Then by recalling 
that the left- and right-handed sectors of the theory approximately decouple, we could see that all the massive diquarks 
come in pairs of degenerate parity-even (scalar) and parity-odd (pseudoscalar) states. 

Regarding the nature of the massive diquark states, let us note that they may truly be just resonances in the full 
theory, since they could decay into the pseudo-NG bosons and/or gluons of the unbroken SU{2)c- At high density, 
however, both the running coupling as(/z) and the effective Yukawa coupling gy = |Aq"|/F'^'' ~ |A^|//z are small, 
and, therefore, these massive resonances are narrow. 

At the end, we would like to add a few words about the higher spin channels that we do not study here. In view of 
studies in Ref. 1 13|, it would be of great interest to investigate also the case of the spin one diquarks, as they might be 
rather light in the color superconducting phase. The general form of the BS equations for such diquarks are exactly 
the same as in Eqs. (m) through (poh. Of course, the structure of the BS wave functions would differ. 
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APPENDIX A: ANGULAR INTEGRATION 



We need to calculate the following traces over the Dirac indices: 



tr 



tr 



7"A('=)7'^A(^')1 = (7'"'(I + ee't) - 2ee'g''"g''°t + ee' 



qp^ +^j^ 



n'^^h^f^r = -g'^'il - ee't) + ( .g^o - e^l) ( g'^" - e'^ 



pM 



|91 



b1 



»"°-^'l)('""-4y 



(Al) 
(A2) 



where e,e' = ±1, t — cos0 is the cosine of the angle between three- vectors g and p, and irrelevant antisymmetric 
terms are denoted by the ellipsis. 

By contracting these traces with the projectors of the magnetic, electric and longitudinal types of gluon modes, we 
arrive at 



Ofi^ tr 



rA(-VAH=2(l-t)^!±^!±^^, 
^ ^ J g^ + p"^ — 2qpt 



Oi'J tr 



7''A(-)7''A(+) =(l + t)- 



q2 _^_ p2 _ 2qpt 

(g - p? + (g4 - Pi? 

' qi j^p2 _ 2qpt + (g4 - ^4)^ ' 



q2 ^p2 _ 2qpt + (g4 - J54)2 ' 



(A3) 



Oil) tr f7^A(-)7^A(+)l = 2(1 + t / !/ -'^^1'^ - (1 + „ Jl;<„ti'!\: "^^!.. , (A4) 



(A5) 
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0« tr [7^A(-)7"^A(- 



O^l^ tr 



7'^A(-)7''A(- 

O^^} tr [7^A(-)7''A(- 

OW tr [7'^70A(-)7''7"A(- 

O^l' tr [7-7"A(-)7''7°A(- 

0(3) tr [7'^70A(-)7''7°A(- 

7^7°A(-)7'^7°A(+ 

7''7°A(-)7'^7°A(+ 

O^^} tr 7^70A(")7"^70A(+ 



0«tr 



O^} tr 



2(1 + t) 
2(1 -t) 
(1-i) 



q2 ^ p2 _ 2qpt 
q^ + p^ + qp{l - t) 



ql ^ jp. __ 2qpt 

{q + p? + (qi - Pi? 



(i-t) 



jq + P? + (qj - Pi? 
q2 j_p2 _ 2qpt _|_ (g^ _ p^)2 ' 



q2 _^_p2 _ 2qpt + (^4 - p4)2 ' 
g^ +p^ + gp(l - t) 



-2(1 -i) 2^2 o , ' 

qz ^ pz _ 2qpf 
2qp{l - e) ^ {q - pf - (94 - Pif - 2i{q4 - Pi){q - p) 



q2 _^_p2 _ 2qpt q"^ +p'^ — 2qpt + {q^ — ^4)^ 

{q - P? - {qi - PiY - 2i(g4 -Pi){q-p) 



-(1 + 0- 



g2 + p2 _ 2qpt + (g4 - ^4)2 
q2 _|_p2 _ gp^]^ _|_^^ 



= -2(l + i) 2^ 2 9 , ' 
(7^ + p^ — Iqpt 

2qp{i-e) ^,^_Aq + p?-{qi-Pi?~Hqi-Pi){q + p) 



q2 j^p2 _ 2qpt q'^ +p^ — 2qpt + {q^ — p^Y 

iq+pf - {qA - PaY ~ 2i{qi - pi){q + p) 



-(1-0- 



q2 ^p2 _ 2qpt + (^4 - p4)2 



where q ^ \q\, p = \p\, q^ ^ -iq^ and P4 = -ipo- 
Therefore, 



dOP^,(g-p)tr[7^A(,-'7''AW 



2iTr 
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(2m)= 



In- 



(2m)' 



.3 |e^|3 + 7rM^uj/2 (e^)2 + 2Af2 + c^2 

/r- = q' I dnV^,{q~p) tr [7^A(")7'^A(^)" 



2iTT 



hn. 



{2^Jif 



-^In- 



(2^)^ 



_3 |e^|3+7rAf2^/2 tt (e-)2 + 2Af2 + ^2 
where Af2 = 2as/i2/7r and u; = |(74 — P4|. 

J^— - g2 y d^V^,{q - p) tr [7'^70A(-)7''7°A(-) 

{2pf 



^In 



(2^)2 



(eg-)2+cc;2 



2i7r 



-^h.- 



hi- 



(2m)^ 



3 |e^|3+7rAf2^/2 (e^)2 + 2Af2+, 



-e 



/2-+ = g2 I dnV,,{q-p) tr 7^7°A(-)7'^70A(+) 



2i7r 



-Im. 



(2/i)= 



In- 



(2m)= 



3 \eq\3+TrRPuj/2 [eqY + 2NP + u' 



-^In-^ 



{2^,f 



{eqY+oj^ 



(A6) 

(A7) 

(A8) 

(A9) 

(AlO) 

(All) 

(A12) 

(A13) 

(A14) 



(A15) 



(A16) 



(A17) 



(A18) 



APPENDIX B: A NON-PERTURBATIVE CORRECTION TO THE SD EQUATION 



In hght of our analysis in Sec. IV, one could argue that the SD equation might get a large non-perturbative 
contribution, coming from the pole contributions in the full vertex function, see Eq. (^2|) and Fig. ^. We remind that 
the pole structure of the vertices is related to the existence of the NG and pseudo-NG bosons in the theory (for more 
on this see Sec. VI). 

If one recalls that the SD equation is quite sensitive to the long range dynamics (|P| <C /Lt), it would be very natural 
to ask whether the pole contributions to the vertex function in Eqs. (p3) could modify the SD equation and its 
solution. The revealed non-perturbative contributions could conveniently be combined in the matrix form as follows: 



25 



dr^^'{q + P,q)\ 



P'^Pu 



35^XiT')J' 



AgV- - A,P^ 



-S^^{T%'A,V- 



K^^) 



A\ 3 



A,n 













^^ =AoP_ 



(Bl) 



It is of great importance to notice that this is longitudinal, i.e., 6r^f^{q + P,q) ^ P^ = (P^, P /?>)■ As a result of this 
property, the contraction of this vertex with the transverse (with respect to P) projector of the magnetic gluon modes 
is equal to zero. Regarding the other two types of the gluon modes (electric and longitudinal), the corresponding 
contractions are non-zero, and they lead to a finite contribution to the right hand side of the SD equation (^l|). We 
stress that the vertex in Eq. (Bl) is longitudinal with respect to P^ (notice tilde), while the projector of the electric 
modes is transverse with respect to P^ (no tilde here). This difference is responsible for a non-zero contraction 
involving the electric gluon modes. It is still, however, the longitudinal gluon mode that, after being contracted 
with the vertex in Eq. (Bl), gives the most significant contribution to the SD equation. By performing the explicit 
calculation, we arrive at the following extra term to the right hand side of the gap equation: 



5A; 



TT 
-— ( 

2 



d'^q 



(90 + e„ )A„ 



{q-py 



V (2^)4 ^^ '^ P> ql-{e^r-\A-\Hq-pnq-p)x 



T^f^i^iq-p), 



(B2) 



which results in the following term to the equations for A 



sa: 



167r2 

as 
16^ 



dq^dq 



dqiA^ 



ql 



lAoP 



< + o 



-^+^1^ 



'ql + {^^? 



In- 



(2m)^ 



qi + i^^y 



(B3) 



This correction is of the same order as the correction from the longitudinal gluon modes in the gap equation ( |29| ) 
when the bare vertices are used [0|. Therefore, such an additional correction could only modify the overall constant 
factor in the solution for A^ . The exponential factor and the overall power of the coupling constant in the solution 
should remain intact. 

Of course, as we discussed earlier, the overall constant in the expression for the gap might get other kinds of 
corrections which have not been analyzed here |15| , |31[] . Sorting out all such corrections is a rather complicated 
problem that is outside the scope of this paper. 



APPENDIX C: THE ANALYTICAL SOLUTIONS TO THE BS EQUATIONS 



In this appendix, we present the approximate analytical solutions to the BS equations. In general, our approach 
here resembles the method commonly used for solving the gap equation pl]-p^ . One of the purposes of the analysis 
below is to illustrate that, while the Meissner effect is irrelevant for the pairing dynamics of light bound states with 



M <^ |Aq I, it is crucial for the pairing of quasiclassical bound states (see Subsec. VII A) 



1. NG bosons 



Let us start from the BS equation ( [73| ) for the (pseudo-) NG bosons in the x-doublet channel. After performing 
the standard approximations that were extensively discussed in many studies of the gap equation |ll[-|l6[| (see also 
Sec. [II), we arrive at the following integral equation: 



xiPi) 



2a. 
9tt 



dqixiqi) 



In- 



A 



ql 



lAp-p 



\pi -(74I +c|An 



(CI) 



where, for brevity of notation, we use x = Xi ^-^d A = {AttY'/'^ ^/aj . Notice that, in accordance with Eq. (EJ), the 
Meissner effect is taken into account by the term c|Aq | in the logarithm. Without loss of generality, we substitute 
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c = 1 in what foUowsy. The integral equation (CI) could be approximately reduced to the following second order 
differential equation: 



{P4+\A^\)X"{P4)+X'{PA} + 



xipd 



0, where v 



I Pi -r 1^0 

subject to the infrared and the ultraviolet boundary conditions, 

X'(0) = and x(A)-|Ao-|x'(A)«0. 



8a. 
97r ' 



(C2) 



(C3) 



In order to get the estimate for the solution, we consider the differential equation on the two adjacent intervals, 
< P4 < |A(^| and |A(^| < P4 < A, separately. The approximate analytical solutions that satisfy the boundary 
conditions in Eq. (|C3|) read 



x{Pi) = Ci cxp 



Pi 



2|An 



V 1 — t^-^ cosh 



Pi 



2|A;7l 



— sinh 



Pi 



vT^^^ 



2|An 



, for 0<p4<|Ao 



x(p4) = C2sin('^ln^^^^ 
2 Pi 



for |Ao|<p4<A, 



(C4a) 
(C4b) 



where Ci (z = 1, 2) are the integration constants. By matching the solutions at p4 = |Aq |, we arrive at the following 
condition that determines the value of the gap: 



In- 



A- 



|A- 



= — arctan 

V 



VI - i^2 coth(Vl - i^V2) - 1 



(C5) 



This leads to the same (up to an overall constant of order one) expression for the gap as in Eq. ( p2D . Recall that the 
latter was derived without taking the Meissner effect into account. We conclude, therefore, that the solution for the 
X-doublet NG boson is not sensitive to the screening due to the Meissner effect. 

The analysis of the BS equation for the (pseudo-) NG bosons in the A-doublet and ?7-singlet channels is very similar 
to the analysis for the x-doublet and we do not repeat it here. 



2. Massive diquarks 



Let us consider the integral equation (117) for the BS wave function of the massive x-doublet. It is instructive to 
start with the analysis of this equation by neglecting the Meissner effect at first. This is achieved by substituting 
c = 0. In this special case, the differential equation reads 



PiXl'{Pi) + xl{Pi) + -jK^^^ {pi)x{Pi) = 0, where v 



97r ' 



(C6) 



along with the same boundary conditions as in Eq. (C3). The kernel [compare with Eq. (118)] is approximately given 
by 



K'^^\pi) = 



r |A-|/ (|A-|2 - M2) , for < P4 < zl\A^\, 



iM, 



for zij\Aa \ < Pi < A, 



(C7) 



where zm = ^\A^\^ - M^/IA^]. 

The analytical solutions to the differential equation in two qualitatively different regions are given by 



''The analysis for c 7^ 1 is a little more complicated since the two scales, c|Aq | and |Aq |, are different. Despite this, the final 
result for c 7^ 1 would remain qualitatively the same as soon as c is a constant of order one. 
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Ml^O I' 



C2 sin ( - In — 
2 P4 



By matching these two solutions at p4 

A 



-M\^0 I' 



for z 
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M 
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'-Ml'-^O 



= — arctan 
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(C9) 



By comparing this with the gap equation (136) (see also Appendix B in Ref. [lT|), we derive the following spectrum 
of massive diquarks in absence of the Meissner effect: 



M„ = lA" 



exp 



27rn\ 



n = l,2,. 



(CIO) 



Below we argue that none of these massive states survive after the Meissner effect is taken into account. In fact , thi s 
is almost obvious when we notice that the BS wave functions that correspond to the states with masses in Eq. (CIC) 
have rather rich node structure in the far infrared region < p ^ \ Aq \ . Topologically, the n-th wave function has 
exactly n zeros. These n zeros appear at 



pi'' - |Ao 



-^f.fc- arctan 44 



fc = l,2,. 



(Cll) 



In the weakly coupled theory, 1/ <C 1, we find that p\ <$^ \Aq \ for any k. This suggests that after taking the Meissner 
effect back into consideration, the mentioned structure of the nodes in the BS wave functions of the massive ^-diquarks 
would become impossible due to strong screening effects of the gluons in the infrared region < p4 < \Aq\. 

To substantiate the claim of the previous paragraph, let us now consider the equation where the Meissner effect is 
qualitatively taken into account. We arrive at the following differential equation: 



(Pi 



\A^\)x"{pa)+x'{P4) + -K^''\pa)x{pa) = 0, where ^^\j^- 



8as 



(C12) 



The wave function should again satisfy the boundary conditions in Eq. (|C3| ). To get the estimate for the solution, 
we consider the differential equation on the three adjacent intervals, < p4 < z^^|A(^|, z^jJAq"! < P4 < |Aq | and 
\^o\ < P4 < A, separately. The corresponding analytical solutions read 



xiPi) ^ Co exp 



Pi 



2|An 



■M 



— v^ cosh 



Pi\f'- 



M 



2zAf|Ao 



—ZM sinh 



Pi\f- 



M 



2zm\A^ 



for Q<pi<zl,\A, 



Ml^O l> 



x{pi) = Ci 






Pi 

|A^ 



1-x 



+C2-^^iFi ( 1 
|Aol 

X[Pi) = C3 sm - In 

V2 Pi 



-,2,- 



Pi 



for 



-Ml'-^O 



<Pi<\A 



h 



for |A^|<p4<A, 



(C13a) 

(C13b) 
(C13c) 



where G^2 i^ the Meijer's G- function. The solutions in the first and the last regions are chosen so that they satisfy 



the boundary conditions in Eq. (C3) 



We note that the ultraviolet asymptote (C13c) of the BS wave function of a massive doublet coincides with that of 
the NG boson (C4q). Moreover, this property is shared by all massive states that exist, irrespective of the value of 
their mass. Now, unlike the wave functions of (pseudo-) NG bosons which have no nodes, the BS wave functions of 
massive excitations should have at least one zero somewhere in the region < P4 < A. Since, in agreement with our 
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previous statement, there cannot be any nodes in the ultraviolet region |A^| < p4 < A, they might occur only in the 
infrared, < p4 < \A^\. 

By matching the logarithmic derivatives of the solutions at p4 = z^^j Aq | and p4 = | Aq" |, we obtain the two different 
expressions for the integration constant C2' 



Ci = T[zm) and Ci = T{\), 
where the explicit form of function T{z) reads 
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(C15) 



Notice that we used the gap equation ( p5D to derive the second expression in Eq. ( pl4 ). 

The spectrum of massive excitations (if any) should be determined by the solutions of the equation J-{z) = ^{1) 
where z < 1. Note that the obvious solution z — I corresponds to the no- node wave function of the NG boson. By 
studying the equation numerically, we find that there are no solutions which would correspond to wave functions with 
nodes. However, in addition to z = 1 solution, there is another solution for z < 1. This latter also corresponds to a 
wave function without nodes in the whole region of momenta < p4 < A. In fact, its shape barely differs from the 
wave function of the NG boson. In the spectral problem at hand, however, one does not expect having two solutions 
with the same no-node topology. Therefore, we believe that the extra solution is an artifact of the approximations 
used. Its appearance apparently results from two different splitting of the whole region of momenta into separate 
intervals for M = and M 7^ cases. This is also supported by the observation that, because of the Meissner effect, 
the BS wave function in the doublet channel is always almost a constant function (and, therefore, cannot have zeros) 



in the infrared region < P4 < |Ag |. This is seen already from the integral version of the BS equation (117). It is 



natural that there is only one no-node solution to the BS equation. Since the BS wave function of the NG boson is 
such a solution, no other non-trivial solutions should exist in the doublet channel. 

The analysis of the BS equation for the A-antidoublet is similar and we do not repeat it here. The analysis for 
the singlet is presented in Subsec. VII B in detail. There the Meissner effect is qualitatively taken into account by 
considering only the interaction that is mediated by the gluons of the unbroken SU{2)c subgroup. 
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FIG. 1. The diagrammatic representation of the SD equation. 
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FIG. 2. The diagrammatic representation of the BS equation. 
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FIG. 3. The definition of the decay constant. 





FIG. 4. The pole contribution to the vertex as P ^ 0. 
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